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Abstract 

We investigate the dynamics of a single semifiexible filament, under the action 

of a compressing force, using numerical simulations and scaling arguments. 

The force is applied along the end to end vector at one extremity of the 

filament, while the other end is held fixed. We find that, unlike in elastic rods 

the filament folds asymmetrically with a folding length which depends only on 

the bending stiffness n and the applied force /. It is shown that this behavior 

can be attributed to the exponentially falling tension profile in the filament. 

While the folding time tq depends on the initial configuration, the distance 

moved by the terminal point of the filament and the length of the fold scales 

as T^/^ at time r >> tq and is independent of the initial configuration. 
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Recent developments in micro manipulation techniques using optical tweezers, micro 
pipette and high resolution optical imaging have made the study of single filament dynamics 
experimentally feasible • Today it is possible to apply forces of the order of pico Newton 
with precision. This has spurred lot of activity in the study of single bio-polymers. These 
techniques have been used extensively in the study of stretching elasticity and force extension 
relation of bio-polymers Some of these techniques can also be applied to study the 

effect of compressing forces on single polymer filaments @J^. 

In a physical description of their statistical properties, the chemical details of polymers 
can be summarized using a set of parameters. These parameters characterize the individual 
features of each polymer species. Two of the most important parameters are persistence 
length and bending stiffness. Persistence length, Ip, is the distance along the contour of 
the filament over which the tangent-tangent correlation is destroyed by thermal fluctua- 
tions. Semifiexible polymers have contour length comparable to its persistence length. This 
means that unlike in the case of flexible filaments or rigid rods the dynamical properties of 
semifiexible filaments are determined by both conformational entropy and bending stiffness. 
Thermal fluctuations also play an important role in determining the viscoelastic behavior of 
semidilute solutions of semifiexible filaments 

From the biological point of view, a study of the mechanical and statistical mechanical 
properties of semifiexible filaments are important to understand the dynamical behavior of 
cytoskeletal network. Cytoskeleton, which controls the shape recovery and dynamic cellular 
reorganization of living cells is a complex network of polymers Q. At the length scale of 
a cell most of the polymers which build up this network are semifiexible 0. Microtubules, 
Ip = 5.2mm and Actin filaments, Ip = 17/im are some of the examples. 



For length scales much smaller than the persistence length the filaments are like rigid 
rods. Under the action of longitudinal compressing forces these filaments exhibit the classical 
Euler buckling ||I2|. The dynamics of this buckling has been studied recently |T^. In the 
case of bio-polymers like microtubules, such buckling is also observed as a result of the 
force generated by polymerization ^ and by the movement, on them, of kinesin motors 



1^ . On the other hand in flexible polymers with persistence length much smaller than the 
contour length, the compression forces applied at the extremities do not propagate along the 
polymer. This naturally raises the question, what happens when we apply a longitudinal 
compressing force on a semiflexible filament? Here we investigate the dynamics of a single 
semiflexible filament under the action of an external compressing force. 

In this paper we look at the folding of a semiflexible fllament under the action of a 
compressing force acting along the end to end vector. Previous calculations of folding in 
semiflexible polymers assumes uniform tension in the fllament ||T5| , [l3| . We show that, at 
flnite temperature, the tension proflle in a semiflexible fllament in equilibrium decreases 
exponentially along the contour. This tension proflle is shown to have non trivial effects 
on the folding dynamics of filaments. The filament is described by the position vector 
r(s), parameterized by the contour length s of the filament, where < s < L. Variety 
of boundary conditions are possible. Throughout this paper we will discuss only the case 
wherein the curvature of the filament at the two extreme ends are zero (hinged boundary). 
This implies free rotation of the tangent vector at the ends. The s = tip of the filament is 
allowed to move only along the direction of the end to end vector which is taken to be along 
the The position of the s = L tip, r(L), is fixed. An external force f is applied 

at the s = tip along the x— direction. A detailed discussion of other possible boundary 



conditions will be published elsewhere |]I6[. The elastic coefficients for the transverse and 
longitudinal fiuctuations of semifiexible filaments are not the same. This implies that, for 
the above boundary conditions, a compressing force could lead to either folding or buckling. 
However, the transverse spring coefficient in the stiff limit is proportional to stiffness k and 
the longitudinal spring coefficient proportional to k^/T Thus an applied force on a 

filament with hinged ends will always lead to folding. 



Actin filaments exhibit a persistence length of about 17/im at room temperature [jTT 
Critical force for buckling is then of the order of a few pico Newton. In a laser tweezers 
experiment we can apply forces of the order of pico Newton with precision. Typically, in 
such experiments the filament is held by a bead of size ~ 0.5/xm to a precision of ~ ±50nm 



p!8| . In view of these, it should be possible to observe folding dynamics in laser tweezers 
experiments using actin filaments of length < lO/xm or microtubules of length > lOO/^m . 



In our numerical simulations, semiflexible filament are modeled as rods of equal length 
hinged together linearly. Position of each hinge is defined by the Cartesian coordinates 
Fj = {xi,yi). The total energy of the filament, with the compressing force f acting at ri, is 

^ N-2 N-1 

F = -J2il~U- t,+i) + 9^t^ - a^iex ■ f , (1) 
^ 1=1 1=1 

where ti is the unit vector pointing from i to i + 1. The first term is the contribution 
from the curvature of the polymer, k being the elastic constant, and the second term is to 
enforce the local inextensibility condition, ti^ = 1, of the rods. We work in the over damped 
hmit where the dynamics of the filament is described by the Langevin equations: 



i^Xi dF 

Here the effect of the solvent enters only through the inverse friction coefficient F. Equa- 
tions ^ are solved along with the constraint = 1. The noise r] obeys the fluctuation 



dissipation theorem which in the case of filaments with equal monomer mass leads to |^ 

Integrating equations ^ we get the increments Axi and Ayi, which are functions of 
the unknown Lagrange multipliers Qi. Increment in time is represented by the dimensionless 
variable At = TAt . The coordinates of the vertices at time r + Ar obtained from the above 
equations have to satisfy the constraint t^(r + Ar) = 1. These A^ — 1 simultaneous nonlinear 
equations are solved to determine the Lagrange multipliers gi. We choose Axn = 0, Ay^ = 
and Ayo = at all times. There is no constraint on the angle between the filament and the 
force direction at both ends. 



We have carried out extensive numerical simulation of the above model using a chain 
with length = 150 for different values of n and /. One run in the simulation correspond to 
At = 1. All results are averaged over 10 initial configurations. The initial conformation of 
the filament (r = 0) is chosen from an equilibrium distribution which satisfy the boundary 



condition yi = and Un = The snapshots from our simulation for a constant applied 
force / is shown in figure 1. We observe that the filament first develop a fold at a distance 
If from the end at which the force is applied. // is the position of the first maximum in 
curvature along the contour when the free end is perpendicular to the end to end vector. 
For a given k, and force /, this "folding length" (Zj) is independent of the initial configuration 
and is found to depend on bending stiffness as // ~ y/n. This is shown in figure 2. 
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FIG. 1 The snapshots of conformations of a semifiexible filament with k = 500 and 
N = 150, taken at equal intervals of time At = 0.2 million for an applied force / = 30. 
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FIG. 2 A plot of log^o \ogiQlf. The filled squares are from the simulation. The 
dotted line is a fit which shows // ~ y/n. 

A better understanding of the dependence of folding length on k is obtained from the 



continuum equations for the worm like chain model pT| , p2| . The equation of motion in the 
over damped limit for the boundary conditions, r(L) =a constant and y{0) = 0, is then 

+ (/ + 2g(s)r(s) - !,;^r) + ,,{s,t). (3) 

here dots denote the derivatives with respect to arc length. ri{s,t) is the thermal noise and 
r is a kinetic coefficient. Using the constraint r(s)^ = 1, and equation |^ we get, 

1 (9f2 

-— = 2g - 2gr' - l,r . d',T = 0. (4) 

Comparing the curvature and force terms in equation ^, we get a length which scales as 
I ~ \J\k/ f). Similar length scale is also obtained from the constraint equation. 

Let us try to get some insight into the tension profile along the filament. Dropping higher 
order derivative terms from the constraint equation ^ leads to 

— — =gf - g^r^ = {). (5) 
2V dt ^ ^ ^ ' 

We use g^ to distinguish this tension from that obtained with k. Numerical solution of 
the above equation for time t = is shown in figure 3. Here we have replaced r(s)^ by its 
approximate ensemble averaged value 

s / 1 \ EnSin^f^l 
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FIG.3 g{s) obtained from equation | for k = 500, 1000, 1500, 2000, 2500 and n=150. 



We see that the resulting tension decays exponentially. We find the length scale to vary 
as a/k. The tension profile is found to depend on the number of modes taken in equation ^ 
and is exponential only in the limit when the number of modes are large. The lower cut off 
in length being set by the segment length. 

We now compare this tension profile with that obtained from the simulations. In the 
simulations, the tension Qi has contributions from three sources, namely, curvature, noise 
and the applied force. Like in the case of continuum equations we can calculate a tension 
g{ by setting the noise and curvature terms to zero. Since the equations of motion are 
nonlinear one may doubt the validity of a comparison between the tension g{ obtained from 
this calculation with g-^{s) obtained from the continuum equations. To dispel this doubt we 
linearize the equations for the tension profile by neglecting terms of the order (Ax,)^ and 
(A|/j)^, in equation = 1. The results obtained from such a linear calculation agree well 
with the full non-linear results (data not shown). However, this linearization also means 
deviation from the constraint = 1. Even for small At these errors in segment length can 
become very large in a short time. Therefore we have to use the full non-linear equations 
for updating the filament configuration. 

In figure 4 we show this tension profile function of the vertex index i. The tension 

g{ goes to zero exponentially as predicted earlier. The distance ^ over which the tension 
falls to half its maximum value is found to increase as ^/k for a fixed value of /. 
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FIG. 4 / / as a function of the contour length i from the point at which force is apphed, 
for K = 300, 1000, 2000, 3000. Inset: log(/t) is plotted against the log^, showing ^ ~ k^/^ 



We now look at how the buckling length If depends on the applied force /. We have seen 
in the previous paragraph that the tension in a semiflexible filament falls exponentially along 
its contour. This exponential tension profile differentiates between semiflexible filaments and 
stiff elastic rods wherein the tension profile is uniform. Thus it is instructive to extend the 



stability calculations for a rod with uniform tension [12] to the present case of exponential 
tension profile of the form exp(— as/i/zt) + 1. The equation of equilibrium is. 



9^y ld_ 

^ds"" 2ds 



ds 



exp(— as/v^) + 1 



0. 



(7) 



These equations are solved numerically with the boundary conditions y\s=L = 0, y\s=o = 
0, y{0) = 0, y{L) = 0. The lowest energy solutions of this equation show that the value 
s = Smaxi of the first maximum in curvature, shifts to smaller values of s with increasing /. 



This is indeed what we see in our simulations of the semiflexible filament. 
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FIG. 5 Folding length (//) as a function of the applied force(/); The dotted line represents 
the buckling point obtained from the numerical solution of equation |^ with a = 10. The 



straight line corresponds to / 



/ 



In figure 5 we compare the buckling length / f obtained from the simulation with Smax from 
the solution of equation ^ a is the only adjustable parameter used to get the agreement 
between the two results. 

We now come to the dynamical behavior of the fold. The length of the segment before 
the fold, the end to end distance and tension g{ all depend on time. As can be seen from the 
configuration snap shots in figure 1, once the filament has folded the applied force acts to 
pull the filament. This implies a change in sign of tension g{ . At r = we have g( = —f ■ 
As the filament folds g( changes sign saturating to g( = f . The time tq at which (^i (tq) = 
can be designated as the folding time. We find this value of the folding time to depend 
on the initial configuration of the filament. The end of the filament at which the force is 
applied, moves by a distance Xi{r) at time r. This is depicted in figure 6 along with the 
length of the fold. 

For T >> tq all the distances scale as r^/^. This shows that in this regime the problem 
is similar to that of stretching a tense string ||2^. For large / we can then neglect the 
contribution from the curvature energy and noise. The dynamics is then determined by the 
straight segment up to the fold moving through the fluid. This results in a linear tension 



profile along the filament |22], with the tension falling to zero approximately at the fold. 
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At s = 0, g^{s = 0) = /, by neglecting the curvature and noise terms on the right hand side 
of the equation ^ using dimensional analysis, we get, x{s = 0, r) ~ (t)^^^- This is shown 
in figure 6. 
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FIG. 6 The displacement of xi{t) after folding of the filament (Filled squares). Crosses 
shows the position of the fold along the contour. Open circles represent the arc length up 
to which g-f is linear, (r — toY^"^ line (solid line) is given as a guide to the eye 

As the end moves, the position of the fold moves along the contour ( see figure 1). Hence 
the fold length and the point at which tension falls to zero too follow similar behavior as 
that of x{s = 0, r). 

To conclude, we have studied the dynamics of a single semifiexible filament under the 
action of a compressing force at one end which lead to folding. We have shown that the 
folding length(/j) of the filament is proportional to ^/K, and decreases monotonically with 
the applied force /. This follows from the fact that the tension in semifiexible filaments 
decay exponentially along the contour. It is shown that buckling of a homogeneous elastic 
rod with an exponential tension profile is qualitatively similar to that of the folding of a 
semifiexible filament. The folding length for a given k, and force f does depend on the 
strength of the noise t]. Thus once the initial configuration is chosen, a "zero temperature" 
calculation should give the same results. As the filament folds the tension changes sign and 
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we move into the stretching regime. We have shown, using numerical simulations and scaling 
arguments, that in this regime the fold length, the end point Xi{t) and the distance along 
the contour to which tension has propagated, scales as r^^^. 
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